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Abstract 

We consider weakly regular Gowdy-symmetric spacetimes on T 3 satisfying the Einstein- 
Euler equations of general relativity, and we solve the initial value problem when the 
initial data set has bounded variation, only, so that the corresponding spacetime may 
contain impulsive gravitational waves as well as shock waves. By analyzing, both, future 
expanding and future contracting spacetimes, we establish the existence of a global foliation 
by spacelike hypersurfaces so that the time function coincides with the area of the surfaces 
of symmetry and asymptotically approaches infinity in the expanding case and zero in the 
contracting case. More precisely, the latter property in the contracting case holds provided 
the mass density does not exceed a certain threshold, which is a natural assumption since 
certain exceptional data with sufficiently large mass density are known to give rise to a 
Cauchy horizon, on which the area function attains a positive value. An earlier result by 
LeFloch and Rendall assumed a different class of weak regularity and did not determine 
the range of the area function in the contracting case. Our method of proof is based on a 
version of the random choice scheme adapted to the Einstein equations for the symmetry 
and regularity class under consideration. We also analyze the Einstein constraint equations 
under weak regularity. 



1. Introduction 

We study here the class of Gowdy-symmetric spacetimes on T 3 when the matter model 
is chosen to be a compressible perfect fluid. We formulate the initial value problem when an 
initial data set is prescribed on a spacelike hypersurface and we search for a corresponding 
solution to the Einstein-Euler equations. In comparison with earlier works on vacuum 
spacetimes or on the coupling with the Vlasov equation of kinetic 



theory JMQIIJ 28], dealing with compressible fluids is comparatively more challenging, 
since shock waves generically form in the fluid \24 \ and one must work with weak solutions 
to the Einstein-Euler equations. 



In the present paper, building upon work by LeFloch and co-authors in [13]-|19l]. we 



deal with the initial value problem associated with the Einstein-Euler equations under 



1 Laboratoire Jacques-Louis Lions & Centre National de la Recherche Scientifique, Universite Pierre et Marie 
Curie (Paris 6), 4 Place Jussieu, 75252 Paris, France. Email: grubic@ann.jussieu.fr, contact@philippelefloch.org. 
Blog: philippelefloch.org. 

Accepted for publication in Archive for Rational Mechanics and Analysis. 



weak regularity conditions, only: the first-order derivatives of the metric coefficients 
and the fluid variables are assumed to have bounded variation, which permits impulsive 
gravitational waves in the geometry and shock waves in the fluid. Analyzing, both, 
future expanding and future contracting spacetimes, we establish the existence of a global 
foliation by spacelike hypersurfaces, so that the time function coincides with the area of 
the surfaces of symmetry and asymptotically approaches +00 in the expanding case and 
in the contracting case. More precisely, the latter holds provided certain exceptional initial 
data are excluded which, instead, give rise to a Cauchy horizon ijzjjlsoll . 

Recall that LeFloch and Rendall [16] established the existence of weak solutions to the 
initial-value problem, but assumed a different class of regularity and did not describe the 
interval of variation of the area function in the contracting case. The method of proof 
introduced in the present paper is based on a generalized version of the random choice 
scheme originally introduced by Glimm |Q] (cf. the textbooks U which we need to 
adapt to the Einstein equations and to the (Gowdy) symmetry class under consideration. 
In order to show the convergence of the proposed scheme, we derive a uniform bound on 
the sup-norm and the total variation of, both, first-order derivatives of the essential metric 
coefficients and the matter variables. Metrics with low regularity were considered earlier 
in J3L @] for the class of radially symmetric spacetimes, which, however, do not permit 
gravitational waves which are of main interest in the present work. 

We thus consider spacetimes (At, g) satisfying the Einstein-Euler equations 

Gj = Tj (1.1) 

and study the initial value problem when an initial data set, describing the initial geometry 
and fluid variables, is prescribed on a spacelike hypersurface with 3-torus topology T 3 ; 
moreover, we assume that the initial data set is Gowdy-symmetric, that is, are invariant 
under the Lie group T 2 with vanishing twist constants. (See Section 2, below, for the 
precise definitions.) The left-hand side of (|l.l|l is the Einstein tensor G a p := R a p - (R/2)g a p 
and describes the geometry of the spacetime, while R a p denotes the Ricci curvature tensor 
associated with g. By convention, all Greek indices take here the values 0, . . . , 3. 
The stress-energy tensor of a perfect compressible fluid reads 

Tj:=(p + p)u a u?+pg/, (1.2) 

where p represents the mass-energy density of the fluid and u its future-oriented, unit 
time-like, velocity field. We also assume that the pressure p depends linearly upon p., i.e. 

p := tfp, (1.3) 

in which the constant k 6 (0, 1) represents the speed of sound in the fluid and does not 
exceed the speed of light, normalized to be unit. 

Two main results are established in the present paper, for expanding spacetimes in 
Section [6] (cf. Theorem 16. Il l and for contracting spacetimes in Section [7] (cf. Theorem 17. fj) . 
The main elements of proof are as follows: 

• Our first task is to formulate the initial value problem for the Einstein-Euler equations 
by solely assuming weak regularity conditions on the initial data set. Specifically, 
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in areal coordinates (£, 9), the fluid variables as well as the first-order derivatives of 
the metric coefficients are assumed to have bounded variation in space, in the sense 
that their (first-order) distributional derivatives are bounded measures. This is so 
except for one metric coefficient which is denoted by a (see Section 2, below) and 
describes the conformal metric of the quotient spacetime under the group action of 
the Gowdy symmetry. The time-derivative a t has bounded variation in time, but no 
such regularity is imposed (nor available) for the spatial derivative ag. Cf. Sections |2] 
and|H below, for further details. 

Next, in Section SI we study special solutions to the Einstein-Euler system. On 



one hand, we analyze the Riemann problem (first studied in |29t 12711 , as far as the 
fluid variables in the flat Minkowski geometry are concerned) associated with a 
homogeneous version of the Einstein-Euler system. On the other hand, we study 
here an ordinary differential system which takes into account the source terms, only, 
after neglecting the propagation-related terms of the Einstein-Euler system. Here, a 
particular definition of the "main geometric variables" is essential in order for certain 
linear interaction estimates to hold and in order to avoid any resonance phenomena, 
due to the possibility of vanishing wave speeds. 

A generalization of the random choice scheme for the Einstein-Euler equations is 
introduced in Section |5l which allows us to establish a continuation criterion. In 
short, solutions to the Einstein-Euler equations are proven to exist as long as the 
mass-energy density variable \i (or, more precisely, the geometric coefficient a) does 
not blow-up. 

Based on the above arguments, we construct solutions to the initial value problem 
and investigate the global geometry associated with initial data sets generating future 
expanding (in Section [6]| and future contracting (in Section [7} spacetimes. We prove 
that, in the expanding case, the areal foliation is defined for all values of the area 
function, while, in the contracting case, we establish the same global result for initial 
data with "sufficiently small" mass density. The latter is a natural assumption since 
certain exceptional data (with sufficiently large mass density) are known to give rise 
to a Cauchy horizon. 



2. Einstein-Euler spacetimes with weak regularity 

2.1. Euler equations 

Our first task is to introduce a definition of weak solutions to the Einstein-Euler equa- 
tions. We focus here on the formulation of the Euler equations, and we refer to [17] for a 
detailed treatment of the (vacuum) Einstein equations. We are given a smooth manifold M 
(described by sufficiently regular transition maps) together with a locally Lipschitz contin- 
uous Lorentzian metric g. We suppose that this spacetime is foliated by compact spacelike 
hypersurfaces "TYf (diffeomorphic to a fixed three-manifold e H), i.e. 

M=\jn, 
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determined as level sets of a locally Lipschitz continuous time function t : Ai — > I C R (I 
being an interval) so that g a ^d a t is a future-oriented timelike vector. The lapse function N 
and the future-directed, unit normal n — (n a ) are then defined as 

N := ( - g^djdpty 1 ' 2 , rf := -Ng a %t. 

These fields are measurable and bounded, only, that is, belong to Lj^ (AI). We assume, in 
addition, that all fields in L^ c (Al) actually admit well-defined values on each hypersurface 
fit, which also belong to L°°(Ht)- We write, in short, N, n a e L™ c (L°°), by specifying first the 
regularity in time and, then, the regularity in space. In practice, these objects will also enjoy 
certain continuity properties in time, as we will state it below. Furthermore, we impose 
that the lapse function is also locally Lipschitz continuous — which implies that the volume 
forms of, both, the Lorentzian metric and the induced Riemannian metric are also locally 
Lipschitz continuous. 

Under the above assumptions, we decompose the Lorentzian metric in the form 

g = -N 2 dt 2 + (3) g, 

in which ^g is the induced metric on the slices of the foliation. In local coordinates (t,^?) 
adapted to the foliation, we write ®g = (g a {,) with a, b - 1, 2, 3. We introduce the Levi- 
Cevita connection V a associated with g, as well as the Levi-Cevita connection V„ associated 
with ®)g, and we define the second fundamental form k — k t of a slice 'Ht by 

k(X,Y) := g(V x Xn) = -g(V x n,Y) 

for any tangent vectors X, Y. Observe that, at this juncture, V a , V fl , and k a \, belong to L^ c (L 00 ), 
only. 

We now turn our attention to the Euler equations by first assuming that the fluid 
variables are locally Lipschitz continuous functions (which is too strong an assumption to 
allow for shock waves and will be weakened in Definition l2.1[ below). Using the projection 
operator h a P := g a P + n a n?, we decompose the stress-energy tensor (|1.2[) by setting 

p:=T a V^, f := -V^np, S a ? := V 6 h*hl, (2.1) 

so that 

T a ? = pn a n$ + /V + fn a + S a K 

In coordinates x = (x a ) = (f,x fl ) adapted to the foliation, we have n° - 1/N, n a = 0, f - 0, 
and S 0a = and, in order to simplify the notation, we define co := (det(^'g)) 1 ^ 2 , so that 
| det(#)| 1/2 = N a). 

The Bianchi identities imply the Euler equations V a T^ = 0, which yield 

f T a B V a tfdV g = (2.2) 
Jm 

for every smooth and compactly supported vector field £. Here, the volume form dV g — 
| det(g)! 1 / 2 dx = Ncodx is locally Lipschitz continuous, while the derivative operator V is 

L i:c' onl y- 
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On one hand, by taking the vector field £ to be normal, that is, £ = cpd t for some smooth 
and compactly supported function cp, we deduce from (|2.2|l that 

0= f T a J t d a( pdV g + f T a B V a d p t cpdV g 
JM F JM 1 

= J N 2 (pn a + f)d a cpcodx- \ Npnpn^arfcpcodx- \ N S a Q V a ^ t cp codx, 

JM JM ' JM P 

thus 

J ^Npd f cp + N 2 fd a cp + (pd f N + N 2 S ab k ab )cp^codx = 0. (2.3) 

On the other hand, by taking £ to be tangential to the slices, that is, having E? - dt(£) - 0, 
we find 



we find 

= f T a s V a E?dV, 

JM 1 



kJJVK 

f N (pn«n^ + fnp + + S") V tt ^ 
= - f pd b N£ h dx+ \ Nk ab fl b co 2 dx + \ Nj p n a V a ^codx 

JM JM JM 

+ f NS a b V a Z b todx. 
JM 

The third term on the right-hand side reads 

I Ncojfsn a V a ^dx= I j p d t ^codx+ \ N j^V a n^ codx 
JM JM JM 

= r j^codx- r Nk ab fz b codx 

JM JM 

and, therefore, 

J M (jbdd h + N S a b V a E, b - p d b N £ 6 j codx = 0. (2.4) 

At this juncture, we observe that d2.3l )-(|2.4 |) make sense under a weak regularity as- 
sumption on the fluid variables, which motivates us to introduce the following definition. 

Definition 2.1. Consider a spacetime (M, g) with locally Lipschitz continuous metric endowed 
with a foliation whose lapse function N is also locally Lipschitz continuous. Let p > be a scalar 
field and u = (u a ) be a future-oriented, unit timelike vector field, both of them having solely 
regularity. Then, (p, u) is called a weak solution to the Euler equations (|1.2|) if the two identities 

f J (N p{p, u) d t cp + N 2 f{p, u)d a cp + [p{p, u) d,N + N 2 S ab (p, u)k ab ) cp) dtdV< 

f f (j b (p, u) dt? + N S a b {p, u) V a E, h - pip, u) d b N E b ) dtdV (3) = 0, 
Jtel J<H '' /a 

hold for all smooth and compactly supported fields cp and E, a . 



f(3)a - 0, 



In other words, observing that dtdV(3) g = co dx, the conditions in the definition mean 
that 

d t (No) p) + d a (N 2 co f) = co(p d t N + N 2 k ab S ab ), 

d t [coj b ) + coV a (NSl) = -copd b N 

hold in the distribution sense in the coordinates (£,x"). This latter statement does make 
sense since the product coV a (N S?) is understood geometrically in a weak sense as precisely 
stated in the above definition. Our assumption that the lapse function N is locally Lipschitz 
continuous implies that the terms d t N and df,N are bounded on each compact time interval 
so that the right-hand sides of 02.61) are well-defined as locally bounded functions. 

From now on, following [16], we restrict attention to spacetimes with Gowdy symmetry 
on T 3 . Introducing coordinates (t,x ,x B ) (with B - 2,3) that are compatible with this 
symmetry so that (x 2 , x 3 ) e T 2 correspond to the 2-surfaces generated by the Killing fields, 
we can now rewrite the Euler equations. As observed in [16], the Einstein's constraint 
equations (stated below) imply the second and third components of the fluid velocity 
vanish, i.e. 

u 2 = u 3 = 0. 

Consequently, the remaining non-vanishing components of the stress-energy tensor read 

p = N 2 (p. + p) (u ) 2 - p, f =N(p + p) u Q u x , 

S 11 = (p + p)(u 1 ) 2 + pg 1 \ S BC = pg BC . 

The term coV a (N Sf) arising in d2.6l) can be simplified, as follows: 

coV a (N S») = d^NcoS]) - Nco l\ a S a b 

= dtiNcoSl) - Nco (r^ S\ + I c lB S B C ), 

where the Christoffel symbols Th, T^ B e L™ c (L°°) are given by 

We can finally write the Euler equations in the form 

d t (Nco p) + d^rfcoj 1 ) = Nco E 0/ 
d t (coj 1 ) + d 1 (Na}Sl) = coL 1 , 



(2.7) 



where 



E :=N(k u (n +p)(u 1 ) 2 +pttk) + (N(p +p)(u ) 2 -p/N)d t N, 
Ei :=-(N 2 (p. + p)(u ) 2 -p)d 1 N+^Npg BC d 1 g BC (2.8) 

+ \Nd 1 g n ((p + p)(u 1 ) 2 + pg n ). 



According to Definition 12.11 the equations (|2.7[) -(|2.8 [) are required to hold in the distribu- 
tional sense in the coordinates (f, x 1 ). 
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2.2. Weak formulation of the Einstein equations in areal coordinates 
We express the spacetime metric in areal coordinates, that is, 

g = e 2 ^~ u) (-a 2 dt 2 + d6 2 ) + e 2U (dx + Ady) 2 + e~ 2U t 2 dy 2 , (2.9) 

which depends upon four scalar functions U, A, rj, a depending on the variables t, 6, with 
a > 0. We assume that 6 e S 1 - [0, 1] (with a periodic boundary condition), and we 
distinguish between two classes of initial data imposed on a slice < Hf and we solve in the 
future direction t > to: 

Future expanding spacetimes: to > 0, 
Future contracting spacetimes: to < 0. 

The global geometry of the spacetimes is markedly different in each case, but yet many 
fondamental estimates are similar (up to a change of sign). At this juncture, we may expect 
that the time variable should describe the interval [to, +oo) in the expanding case and the 
interval [to, 0) in the contracting case, although we will actually see that dealing with the 
contracting case is more involved. 

After a tedious but straightforward computation from the Einstein equations, one 
obtains the evolution equations 

4(1 

(tfl -1 ir t ) t - (t<*Ue) e = e —(A 2 - a 2 A 2 e ) + taU u , 
(r 1 0~%) f - (r 1 aA e ) Q = ^rV^UfAf - a 2 U e A e ) + all A , (2.11) 

{ta- x ri t ) t - (ta(rj + loga) e ) g = It all 2 + C —A 2 + taW(t), 

understood in the sense of distributions, where the lower-order matter terms are given by 

:= I e 2( J? -u) (p _ Sl + S2 _ S3)/ 

IT 4 :=2e 2 ^S 23 , IP := e 2 ^( p - S 3 ), 

in terms of the spatial part S a ^ of the energy-momentum tensor T a p. Here, relevant com- 
ponents of S a p are given with respect to the orthonormal frame 



e o. = a - le - {ri -u )du £ i = e -(n-u) de> 
e 2 :=e~ u d x , e 3 := y(-A9 x + d y ), 

and read 

Si := S(e\ e 1 ) = e 2 ^ (/i + V ){u 1 ) 2 + p, S 23 := S(e 2 , e 3 ) = 0, 
S 2 -S(e 2 ,e 2 ) = p, S 3 :=S(eV) = p. 

By introducing the scalar velocity v := ^ and using (au ) 2 - (u 1 ) 2 = e~ 2 ^~ u \ we can express 
u°, u 1 in terms of v. 

e 20 ? -U) (u )2 = 1 e 2(,-U) (M l )2 - 



(2.12) 



a 2 (l-v 2 )' 1 v ' 1-v 2 ' 
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Therefore, the matter part of the evolution equations simplifies and we have 

Finally, we can state Einstein's constraint equations 

j = (U 2 + a 2 U 2 ) + "—(A 2 + a 2 A 2 g ) + e 2 ^a 2 p, 

1 f = - U f a + 2U t U 9 + e ^A t A e - e^aj\ (2-13) 

^ = _ fl 3 e 2(,,-U) (p _ Si) _ 

Remark 2.2. The "classical" formulation of the Einstein equations in ureal coordinates involves 
the operator d n - a 2 dse applied to U,A, and rj. However, as pointed out in ill 617 . this formulation of 
the matter equations can not be used to handle weak solutions, since the metric coefficients do not 
have enough regularity in the present context and, for instance, the product a 2 Ugg does not make 
sense. 

2.3. Weak formulation of the Euler equations in areal coordinates 

For the derivation of the Euler equations in areal coordinates, we need the following 
expressions of the metric coefficients: 

N = ae^~ u \ g 00 = -N 2 = -a 2 e 2 ^~ u \ 

gn = e^~ U \ g n = e- 2( "" U) , 

det(( 2 >£) = t 2 , co 2 = det(( 3 >g) = e 2 ^ u h 2 , 

tr/c = -fl"V ( ''- u) (y + rjt - U t ), k n = -a- l e^- u {^ - U t ). 
We can then rewrite the right-hand sides 02.81 ) of the Euler equations and obtain 
Lq := Nco E = a)N 2 (k n (p + pXu 1 ) 2 + (trk)p) + co(N 2 (p + p)(u ) 2 - p)d t N 
= - ate^- U) (n t - U t )(p + p)(w 1 ) 2 - ate 2 ^- u \- t + rj t ~ U t )p 

+ ate 2 ^~ u H^ + m ~ U t ){a 2 e 2 ^ u \p + p)(u ) 2 - p) 
= a feZOHi^ _ l v + ( 7?f _ Ut )(p -p) + a j((p + p)^-^ ~ V)\ 



and 



E; := coZ, = te'i~ u ( - (N 2 (p + p)(u ) 2 - p)d e N + ^Npg BC d e g BC 

+ \Nd e g 11 ((p + p)(u 1 ) 2 +pg 11 )) 
= a te 2 ^~ u \ - (ne - U e )(p -p)- ~~((/^ + " v)\ 



8 



Finally, we obtain 

d{a + 0, + P)^)) + d e (a 2 fc«i-% + p)^) = lf , 

d t (te 2 ^ + p)^) + <9 (« te^fe + Ox + p)^)) = E;. 

Recalling the linear equation of state assumed for p, that is, p = k 2 \i, the weak form of Euler 
equations in areal coordinates takes the following final form 

d t (a te^-u)^±J^pj + de ^i ^iHi) M <l±^E J = a te^pL'J, (2.14) 

dU^^^\ + d e [a te^p.t^pj = a te 2 ^^, (2.15) 
understood in the sense of distributions, where Eg , E" are given by 

w/ fc 2 / ,rw. a t l + k 2 v 2 

u t a 1 -v 2 

3. First-order formulation of the Einstein-Euler equations 

3.2. Spacetimes with bounded variation 

For the analysis in the present paper, it is convenient to refomulate the Einstein equa- 
tions as a first-order hyperbolic system supplemented with constraint equations. In fact, 
we will first drop the constraint equations and study the evolution part of this system. 
We emphasize that the Einstein constraint equations are used in the first-order evolution 
equations. 

We now arrive at a new formulation of the field equations as a system of first-order, 
hyperbolic balance laws supplemented with ordinary differential equations. To this end, 
we introduce the new variables 

W± := (U t + aUg), V ± := e 2U (A t + aAg), (3.1) 

which we refer to as the first-order geometric variables, to which we append the fluid 
variables (Jl,v), and the geometric coefficients (fl,v). Here Jl is obtained by scaling the 
energy density and we also introduce a new metric coefficient, as follows: 

v = r\ + log(fl), y. := e 2(v " u) ^. 

Observe that the product containing ag on the right-hand side of the "second" Euler 
equation (|2.15l) may not be defined under the regularity assumptions we are interested in. 
However, if it is multiplied by a, the second Euler equation can be rewritten so that ag 
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cancels out. In this case, using our new notation, the source-terms Eg", E"' of the Euler 
equations take the form 

= E" = --— - -(W+ + W_) + -Ml + W 2 ) + —{V 2 + y 2 ), 

(l-k 2 )t 2 V 7 2 V + 8f + " 

w„ „ v „ ^ «f (1 + k 2 )v k 2 +v 2 

Lj = aE + =- 

fl 1 - w z 1 - u z (3 2) 

= - (v 9 -U e )(l-^) + ^<l^ 

fl \-v l 

= - w_) + ^(w 2 - w 2 ) + i(v 2 - y 2 ). 

Without specifying yet (see below) the regularity of the spacetime metric and lapse 
function, we now collect, in the new notation, the first two evolution equations, the ordinary 
differential equation for a, and Euler equations. 

Definition 3.1. The first-order formulation of the Einstein-Euler equations for Gowdy- 
symmetric spacetimes consists of five evolution equations for the geometry 

1 (W+ - W-) 1 (a 9 

(a-HW ± ) t ± (tW ± = ± + =-7 V + V_ + £ (1 - fc 2 ), (3.3) 

2a 2at 2a 

(fl- 1 r 1 y ± ) f ± {t- x v ± ) e = + {V+ ~J~ ] - | w ±^- ( 3 - 4 ) 

a t = - atji(l - k 2 ), (3.5) 

two equations for the fluid 



2\ v'" 
o ' 



(3.6) 



and f/zree constraint equations: 

v t = V 2 + w 2 ) + i(y 2 + y 2 ) + 

L or 1 — v 

t , 1 9 , _ (1 + /c 2 )o 

-«v e = -(w 2 - w 2 ) + -(y 2 - y 2 ) + ru K - — £, 

2 or 1 — v A 

and 

(ta-\v t + tji{\ - k 2 ))) t - (tave)e =^(W + - W_) 2 + ^(V + + y_) 2 

i ~ (1 + ^ 2 ) 
l-v 1 



(3.7) 



(3.8) 



As we will show it, (|3.7|) are indeed constraint equations, that is, provided the initial 
data satisfy these equations, they then hold throughout the spacetime. This is a standard 
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property; yet, in the present class of weak regularity, it is essential to check it, as we do 
here. 

We are now in a position to state the regularity of the spacetimes under consideration 
and to define our concept of weak solutions. We denote by BV(S 1 ) the space of functions 
with bounded variation, that is, bounded functions / : S 1 — > K whose derivative is a 
bounded measure. The BV regularity for tensor fields is defined by considering each 
component in adapted coordinates. 

Definition 3.2. A foliated spacetime with bounded variation (or BV foliated spacetime, in 
short) is a foliated spacetime (At, g) with locally Lipschitz continuous metric such that the ureal 
metric coefficients Ut, Ue,At,Ae,a,vt,VQ belong to L?° (BV^S 1 )) n Lipi oc (L^(S x )). Such a spacetime 
is said to be an Einstein-Euler spacetime with bounded variation if, in addition, there exist a 
scalar field p and a unit vector field u = (u a ) inE^ c (BV(S 1 ))r\Lipi 0C (L 1 (S 1 )) such that the first-order 
formulation of the Einstein-Euler equations (cf. Definition ^, li holds in the sense of distributions. 

Whenever necessary, we will tacitly extend functions defined on the circle S 1 = [0, 1] to 
periodic functions defined on the real line R with period 1. 

3.2. From first-order to second-order variables 

In this section, we show that it is actually sufficient to solve the set of essential equations 
(|3.3 [H|3.6[) , without imposing the constraints at the initial time. This is justified, as long as 
we can show that there exists a function v satisfying the constraint equations d3.7M3.8D and 
we can recover the metric coefficients U and A from a given BV solution ( W+, V+, a, p, v) to 
the first-order system 03.31 ) - (I3.6D . 

Indeed, by a direct calculation, it follows from Einstein-Euler equations that the right- 
hand side of the constraint equations 03.71 ), that is, 

Cl := L {W 2 + W 2 ) + Kyi + y2_ ) + f ~ + 



2 V + -' 8C + '-' r 1-172' 

C 2 := Uwl - m + I(V* - Vl) + tp^^, 
1 or 1 — v L 

satisfies the compatibility condition 

(fl" 1 C 2 ) f - (Ci) = 0. 

Therefore, the system 

v f = Ci, av e = C 2 , 

can be solved for v, by integrating the first (say) of the above equations 



v(f, 6) := v o (0) + JT ||(Wj + W 2 _) + ^1 + V 2 _) + tp^^\t', 9) dt' . (3.9) 



Obviously, the distributional derivatives v f ,vg of v belong to BV(S 1 ), while v is periodic, 
since vo is periodic and, by construction, 

d t f v e d9= f d t (a- 1 C 2 )d6= \ d e C 1 dd = 0. 
Js 1 Js 1 Js 1 
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Moreover, by a straightforward calculation, it follows that (13.81) holds as well. Observe that 
the above calculations only require multiplying the evolution equations for W+ and V±, 
which is allowed even for weak solutions, since the principal part of these equations is a 
linear operator and jump relations are unchanged. 

It remains to recover the metric coefficients U,A from a BV solution in the variables 
W± f Indeed, by subtracting the evolution equations for W+, we obtain 

(fl _1 (-W + + W_)) + ((W+ + W_)) = 

and, therefore, the system 

U t = ^(W + + W_), U e = Ya { ~ ]N+ + W " } ' 
can be solved in U, by setting 

U{t, 9) = U (6) + \ f (W+ + W_)(f , 9) dt'. (3.10) 
2 Jo 

Clearly, first-order derivatives of U are functions of bounded variation in 9, and U is 
periodic. The same properties hold for the coefficient A, and our argument is completed. 



4. Special solutions to the Einstein-Euler equations 

4.1. Homogeneous solutions 

In this section, we neglect the source terms in the equations d3.3D - d3.6D and we tem- 
porarily assume t = 1, that is tW+ — > W+, t~ l V+ — > V+, tjl — > J[, to facilitate the exposition. 
The results for the original variables are easily recovered after scaling by a suitable power 
of f . Hence, we study the system 

(a- 1 W ± ) t ±(W ± )o = 0, 

(a- 1 V ± ) t ±(y ± ) = O / 
at = 0, 

/ _^ l + k 2 v 2 \ , /- (l + fc> \ n (4-1) 

d *\ a ^-r—r) + M^ 1 _ V 2 1 = ' 



The geometry variables (V ± , W±, a) satisfy linear hyperbolic equations, associated with the 
speeds -a, 0, +a only, whereas the last two equations for ( /I, v) correspond to the system of 
special relativistic fluid dynamics. The fluid equations are coupled to the geometry only 
through a and, thus, can be considered separately from (V+, W±). 

We are interested first in solving the Riemann problem for the above system expressed 
in terms of u :- (W+, V+, a, /i, v), that is, in solving the initial value problem with piece wise 
constant initial data 

u(t ,9)J U! ' l <e °' (4.2) 
u r , 9 > 9 . 
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Here, the constant vectors U\, u r are prescribed within the physical region 



a > 0, y. > 0, \v\ < 1, (4.3) 

and an arbitrary point (tQ, 6q) has been fixed. 

To solve the Riemann problem for the subsystem (a,p,v), we will make use of the 
Riemann invariants y, z, w 

1. (l + v\ k _ 

y : = lo s fl - z:= 2 los l— J + rT^ los ^ 

i k - 

w:= o lo S T-7 -TT7a lo S ^' 



2 1 + ic 2 

This system is not strictly hyperbolic at hypersurfaces v - ±k, where, in each case, two of 
the corresponding characteristic speeds 

v—k v + k 

Aq = 0, \\ = a- -, A2 — a- 



1 — vk' 1 + vk' 

coincide. On the other hand, both Ai and A2 are genuinely nonlinear and \q is obviously 
linearly degenerate. Moreover, we have Aj < A2 where both are strictly monotonically in- 
creasing with respect to v. Up to an additional vector a, they correspond to the characteristic 
speeds of the special relativistic Euler equations, which we denote by 

-~- v + k 

Al,2 



1 + vk' 



Proposition 4.1 (Riemann problem for the Einstein-Euler equations). The system (|4.1|) ad- 
mits a unique self-similar weak solution to the initial value problem (|4.2|) . All the variables contain 
a stationary discontinuity. In addition, W+, V+ contain jumps associated with the speed a r , while 
W-, V- contain jumps associated with the speed -a\. The fluid variables ( Ji, v) have up to three ad- 
ditional waves comprising rarefaction and/or shock waves satisfying Lax shock condition. Moreover, 
the regions 

A M := \{w, z) I - M<w <z< M), M > 0, 

are invariant for the solutions of the Riemann problem associated with the fluid equations. In other 
words, if the Riemann data belong to a set Am, then so does the corresponding solution of the 
Riemann problem. 

The proof of this proposition is given shortly below. We will mesure the strength 
\e{ui,u r )\ of the waves present in the solution of the Riemann problem !/?(w/,w r ) associated 
with Riemann data U\, u r in terms of the strength vector e(u.\, u r ) e R 7 . Its first five compo- 
nents contain the jumps (a, W+, V+), whereas the last two contain the jumps of log Ji across 
shock/rarefaction waves of the fluid equations. The total strength is defined as the sum of 
the components of e, that is 

\e(u h u r )\ := ^ \ei(ui,u r )\. 
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Recall that W+, V+ and fi are all constant across a stationary discontinuity and that a possi- 
ble jump in Jl occurring when one of the rarefaction waves crosses the plane of non-strict 
hyperbolicity would not contribute if included in the definition of e, since p is strictly 
monotone along both rarefaction curves. For the geometry variables, the following state- 
ment is just the triangle inequality, whereas for p the claim follows from the corresponding 
estimate for special relativistic fluids j27lj . 

Proposition 4.2 (Interaction estimates). For any three Riemann problems 91(111, u m ), < R{u m , u r ), 
^(i//, u r ) associated with constant states u\, u m , u r , one has 

\e(ui,u r )\ < \e(ui,u m )\ + \e(u m ,u r )\. 

Lemma 4.3. The strength vector e : Q x Q — > R 7 is a smooth function of its arguments, and for 
all uj, u r , u'j, u' r in a compact subset of CI one has the uniform estimate 



Proof of Proposition \4.1[ We first study the fluid part of the above equations, that is we 
study the subsystem corresponding to u :- (a,Ji,v). Since Ao is linearly degenerate, the 
0-wave curve starting at U\ can be obtained from Rankine-Hugoniot relations or directly as 
an integral curve to the eigenvector Rq := (1, 0, 0). Immediately we have 



The 1- and the 2-wave shock curves S\ t 2 can be obtained by a direct calculation from the 
Rankine-Hugoniot relations, following the exposition in |27l|. Then, it follows immediately 
that a remains constant along both shock curves, but also that p and V are given by the 
corresponding special relativistic shock curves; hence, all of the standard results still hold 
true. For completness, we state S\ respectively S2 starting at U\ in terms of the Riemann 
invariants y, w, z. We have 



e(u\,u\) = e(u h u r ) + 0(l)|e(w/,w,)|(|i/ ; ' - U\\ + \u' r - u r \) 
+ 0(1)IK - «,') - - 



y- = hi, 



V - V\. 



y(u) - y(ui) = 0, 
w(u)-w(ui) = ~ log f + (2Kf3) 





and 



y(u) - y(ui) = 0, 
w(u)-w(ui) = ~ log f + (2Kp) 




2 lo S/-(^ 




where /±(j3) := 1 + jS(l + ^1 + 2//?) is defined for all p > 0, while K := 2k 2 /(l + k 2 ) 2 i 



constant. 



14 



For the shock speeds s,, it can be checked that Si{a,v;a\,V\) - asi(v;vi), where Ti(v;vi) 
have the same form as the shock speeds of the special relativistic equations starting at V\. 
Therefore, Lax shock conditions 

Ai(a,v) < Si(a,v;ai,vi) < Ai(ai,vi), 

immediately follow from the corresponding standard result. On the other hand, for the 
rarefaction curves ^1,^2 starting at U\ we have 

7?! = |u / y(u) - y{u{), z(u) - z{u\), w(u) > w(w/)|, 

%i-\u I y(u) - y{u{), w(u) = w(ui), z(u) > z(u;)}, 

which coincide with the rarefaction curves of the special relativistic equations. Finally, 
the composite wave curves, defined as 'Wi = Si U Ki, are of class C 2 and convex; the 
parametrized curve 'Wi is monotonically decreasing, while IV2 is monotonically increas- 
ing. From now on, we denote the (w, z)-components of IV, by 'Wi. 

The form of the curves 'Wj suggests that they can be regarded as wave curves of the 
special relativistic fluid dynamics; hence, for given Riemann data «/, u r , there exists a 
unique point 

(Jl m ,v m ) e < Wx(]li,vi) n *W z {Ji r ,v r ), 

such that the special relativistic speed of the first wave at v m remains strictly smaller than the 
special relativistic speed of the second one at v r . Therefore, depending on their respective 
signs, it only remains to decide where to put the stationary discontinuity in order to retain 
the correct wave speed relations. In terms of 'Wi, this amounts to jumping from a point 
( \i, v) on the plane a = a\ to the same point on the plane a - a r , whereby all the relevant 
speeds'^, A; prior to the jump get multiplied by a x and all following it by a r . 

For instance, if we have {Jim,Vm) £ S\{Jli, Vi) followed by {Jlr,v r ) e S2{Jim,Vm) and 
< Ti(v m ) < S2(cv) we have to put the stationary discontinuity afirst, i.e., we have four 
constant states 

(ai,Jli,vi), (a r ,Jii,vi), (a r ,Ji m ,v m ), (a r ,Ji r ,v r ). 

Observe that, in the limiting caseTi(v m ) = 0, the middle state (a r , v{) cancels out and the 
solution consists of the three constant states 

(ahjiovi), {a r ,Ji m ,v m ), (a r/ Ji r ,v r ). 

We obtain all the other cases similarly, except if either <R\ crosses the plane v - k or ^2 
crosses v - —k, i.e. if one of the characteristic speeds changes sign. Observe that it is 
not possible for both cases to occur simultaneously, since v > v\ in both cases. If one of 
them occurs however, the solution consists of five constant states, where the stationary 
discontinuity has to be put at v = ±k. If we take, for example, the case of two rarefaction 
waves, where 'Ri crosses the line v - k at some Ji^, then we have 

(fl/,^,y/), (ai,jL k ,k), (a r ,Jlk,k), {a r ,Jlm,v m ), (a r ,JI r/ v r ). 
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We now turn to the geometric variables. The component W+ of the Riemann problem 
consists of at most three constant states where a stationary discontinuity, across which W+ 
is constant, is followed by a jump discontinuity of speed a r across which a is constant. 
We handle W_ similarly and only note that here a jump discontinuity with speed —a is 
followed by a stationary discontinuity. An analogous conclusion holds for V ± and the 
proof is completed. 

Proof of Lemma B31 Since e(ii\, u r ) = if u\ = u r we have 

f 1 de 

e{u\, u r ) = I - — (u/, (1 - t)«/ + zu r )(u r - ui)dz 
Jo our 

and analogously for e{u' ir u' r ); hence, we find 
ei{u\,u' r ) -€i(ui, u r ) 



I 



1 1 de de \ 

I Y~ {Up (1 - T)U\ + TU r ) - ■j — iUl, (1 - T)Ul + TU r ) I (u r - Ui) dz 

- 1 de 



(4.4) 



— (U'j, (1 - T)U\ + XU' r ) (« - U\) - {U r - U,))rfT 
= 0(1) (\u\ - Ui\ + \u' r - U r \) \Ur ~ U,\ + 0(1) \(u' r - u\) - (ll r - U/)|. 

On the other hand, we have \e{u\, u r )\ — 0(l)|w r - ui\, and the claim thus follows. □ 

4.2. Spatially independent solutions 

We set u = (W+, V ± ,a,~p,v) and study the following system of ordinary differential 
equations (ODE, in short) 

u'(t) = g(u(t),t), u(t ) = u 0j (4.5) 

where g represents the right-hand side of the equations d3.3D - d3.6l > and we assume ao, jdo > 
and \vq\ < 1. More precisely, for the geometry variables we have 



1 W+-W 1 ta , 

11 V+-V 2 



(4.6) 



while, for the fluid variables, 

+ ^(W* + + + v^)), 
U^Ci+^Ej' = fl -i f ~ (1 _ ^ _ i (W+ _ w ) 

+ 1 (W 2_ W 2 )+ l (y 2_ y 2 } ) 



(4.7) 



2 V + 8f 
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Finally, the component a satisfies 

a' = -a$f(l - fc 2 ). 
Observe that the above system admits an energy functional 

1 /TA72 , ta /2 \ , 1 nrl , t,2x , 



2 



E := — (W; + Wl) + — x(Yi + VI) + -■ 



2fl v + 8af 2 ' + _/ fl 1-w 2 ' 

whose time-derivative reads 

E > = -I(l ( w + + w.f + -V + - v.f + l {l + k2) 



t\2a ' Sat 1 a 1-v 2 

Hence, £ is non-increasing in (areal) time in the expanding case, while it is non-decreasing 
in the contracting case. We now establish (for the simplified model under consideration in 
the present section, at least) that, as long as the coefficient a remains bounded, the other 
geometric and fluid variables can not blow-up. 

Lemma 4.4 (No blow-up in finite time. Preliminary version). Fix a time to satisfying to > 
in the expanding case and to < in the contracting case. Fix also some initial data uq satisfying the 
physical constraints 

< «o, < Jl , N < h (4.8) 

and consider the corresponding solution of the ordinary differential system d4.5l > defined on its 
maximal interval of existence [to,T). (Here, < to < T in the expanding case and to < T < in 
the contracting case.) Then, the following results hold: 

• In the expanding case, the solution is defined for all time up toT = +oo and 

lim W±(t) = lim = 0, (4.9) 

t— >+oo t— >+oo f 

1 1 

— < ~~7Tr (4-10) 

ao fl(f) 

and 

-l<u<l, 0</.(< +oo. (4.11) 

• In f/ze contracting case, T may be finite; however, as long as a(t) remains bounded, the variables 
W+, V+ remain bounded and fi, v remain strictly inside the physical domain; in other words, 
ifTe (to, 0) and a is defined on the interval [to, T) and remains bounded on that interval, 
i.e. provided 

flo < sup a(t) < +oo, (4.12) 

then the solution u - u(t) extends by continuity to the closed interval [to, T] and, therefore, 
(|4.9|) and (|4.11|l hold in the interval [to, T]. 



In particular, this proposition shows that ( p, v) remain in the physically admissible 
region, i.e. 

< Jl, \v\ < 1. 
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Proof. We will establish here a slightly stronger result than the one stated in the proposition 
and will be able to control the time-dependence of our upper bounds, which will turn out 
to be useful later on. We first rewrite the two fluid equations in terms of log J[ and log j^: 

( log(f £))' = ^0 f(W ± , V ± , v, - (1 - k 2 ) t% 
, l + v\' 2(1 -k 2 ) /TA7 T7 

with 

/(w ±/ v ± , t) = - jt^S - M - ^(w + - w_) + ^(w* - w*) + i(vj - 

+ (1 + v 2 )( - i(W + + W_) + ^(W 2 + W 2 ) + + y2^ 

and 

1 2 1 / -1 

g(W ±/ V ± , v, t) = 777^ - K - 2^ W+ + W ") + 2 (W + + W " } + 87 + 

+ - \w + - w-) + {(wl - w 2 ) + l ( v* _ ^)). 

For the expanding case fo > 0, the function a is non-increasing and thus bounded by the 
initial value. Assume that on some time interval [fo, T\) we have < ju and \v\ < 1. Then, E 
is decreasing and we have that a^W 2 and a~ x t~ 2 V 2 ± are bounded and therefore also W ± , V+. 
Consequently, log Jl and log ^ remain bounded as well, since all the relevant geometry 
expressions are bounded. Therefore, there are constants C,M > depending on T\ and the 
supremum bounds on the geometry such that 

1 _ 

< - < jj < C, \v\ < M < 1. 

Finally, the differential equation for a implies that a is decreasing and remains bounded 
away from zero 

< K < a{t) < a 

for some constant K > and for all t <T\. 

Next, we turn our attention to the contracting case to < and we assume sup te j fo T j |«| < 
+00 for some T < 0. Suppose, moreover, that on the interval [to, T) one has < Jl and \v\ < 1. 
The energy is now increasing and we have the inequality 

E(ty < --E{t), 



which, by Gronwall lemma, implies 



£(0 < E{to)f. 
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Hence, we have 



Wi(t) < a(t)2E(t )-%, Vi(t) < 8E(t )t z a(t), 



'fi' 

that is, the geometry variables remain bounded. On the other hand, this implies as in the 
expanding case that v is bounded away from +1, and p, from and +00. □ 

The bounds in Lemma 14.41 can be stated in a quantitative form, which will be useful 
when dealing with approximate solutions later. We introduce the energy 

(iW+-l/2) 2 (fW_ -1/2) 2 V 2 V 2 _ 



and the energy flux 



h\ := z 1 z 1 + 

2at 2 2at 2 8at 2 8at 2 



(JW+-1/2) 2 (£W_-l/2) 2 V 2 + V 2 _ 



gi : = — 



2at 2 2at 2 8at 2 8at 2 ' 

and observe that they satisfy (for spatially independent solutions) 

at 2ki 
\ = -h - — , 

1 a t 



at gi 
gi = ' 

in which 



(4.13) 



(t(w + + w_) - 1) 2 (v + + y_) 2 

/Ci '.— 1 — 

4at 2 16at 2 

is non-negative, with < k\ <h\. We then introduce the variable^] 

H ± :=a{h x ±g l ). (4.14) 

Proposition 4.5 (No blow-up in finite time). Fix a time to satisfying to > in the expanding 
case and to < in the contracting case, together with initial data Uq satisfying 

< ao, < Jio, \v \ < 1. (4.15) 

Consider the corresponding solution of the ordinary differential system d4.5|) defined on its maximal 
interval of existence [to,T). Then, the following results hold: 

• In the expanding case, T = +00 and for all T > to there exists a constant Cq > depending 
on the times to, V and the initial data uq so that for all t £ [to, T'] 

< H*(t) < max(H-(to),H + (to)), 
a(t) < a{to), 

and 

e -C (t-t ) < < g C (t-fo) 

e -C (f-fo) | z(fo) | < | z(f) | < e Co(t-fo) | z (f )|_ 



1 These variables are natural also for the homogeneous Riemann problem studied in Section |4~T1 and are 
easily checked to satisfy the invariant domain principle, like the Riemann invariants w, z for the fluid; cf . Propo- 
sitionl4~Tl 
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• In the contracting case, T < may be non zero (so that a(t) may blow-up "before" t would 
reach 0), and one has 

tl a{tf , \ 
< H±(f) < - max(H-(t ),H + (t )), 

a(t ) < a(t), 

and w,zalso satisfy the bounds above for all T e (fo, T). 

Of course, controling the two variables h\,gi is sufficient to control the four variables 
W±, V+. Observe also that, in the expanding case, our bounds on H*, w, z are independent 
of any lower bound on a; such a lower bound can however be deduced from the equation 
(loga)f = -(1 -k 2 )t u (and the control of the density /i implied by the above proposition); 
we do not state this here, since such a lower bound can not be "propagated" in the Glimm 
scheme, as only Riemann invariants must be used. 

f 2 a(t) 2 

Observe that, in the contracting case, the coefficient t i a , lo si is greater than 1 and a possible 
amplification of H* may take place. 

Proof. Consider first the expanding case. Here, we have t > and a t < and, since hi and 
ki are both positive, we easily find < h\(t) < /zi(to)- Next, we consider the equation for gi 
and obtain directly (log < so that \g\(t)\ < \g(to)\. After summation, we deduce that 

< H±(f) < —jj-t max(H + (i ),H + (f )), 
<max(H + (f ),H + (f )), 

since a is decreasing. The Riemann invariants w, z for the fluid satisfy a first-order differ- 
ential system whose right-hand side is immediatly controled by H* and, therefore, w, z can 
grow by a multiplicative factor e c ^~ to \ at most. 

Consider next the case of contracting spacetimes. The function is now increasing in 
time and we have 

(fV^iV = -(ht-h)<0, 
v 'a 

since now t < and h\ — k\ > 0. Similarly, one checks that (log (^ 2 fl _1 |^il)) is negative and, 
in turn, we obtain 

i 2 fl(0"^i(f) < %a{t Q )- l h x {to), 

P-a(t)- l \ gl {t)\ < f 2 «(fo) _1 l^i(fo)l, 

which leads to the announced result for H* . The argument for the fluid variables is identical 
to the one in the expanding case. □ 

5. A continuation criterion for the Einstein-Euler equations 

5.1. Generalized random choice scheme 

In this section, we follow [2] and generalize the so-called random-choice scheme to the 
Einstein-Euler system (13.3lHl3.6l ). The scheme is constructed from an approximate solver 
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to the generalized Riemann problem, i.e. the Cauchy problem associated with piecewise 
constant initial data 

u(h,-) = < _ _ (5.1) 

\U r , U > Vq, 

which we denote by < Rg{ui, u r ; to, 60). The approximate solver of interest is constructed by 
evolving the solution u of the corresponding classical Riemann problem < R(ii\, u r ; to, 60) by 
the ODE system (14.6lH|4.7[) . More precisely, we have 

u{t, 9; u,, u r , t Q , Qq) := u{t, 9; u h u r , t , ) + I g(r, S T u(t, 9; u h u r , t , 9 ))dz, (5.2) 

J to 

where S T u"(f, 6; to, 80, u\, u r )) is the solution of the ODE system d4.6l >- (|4.7|l corresponding to 
the initial value u(to) = u(t, 6; U\, u r , to, Oo). 

To formulate the scheme of interest, we denote by r, s > the space and time mesh- 
lengths, respectively, and by (t^, 0^) (for k = 0,1,... and h any integer) the mesh points of 
the grid, that is, 

fjt := to + ks, 9h := hr. 

We also set 

6k,h :=a k + hr, 

where (fljt)jtetsr ^ s a nxe d equidistributed sequence in the interval (-1, 1). We will let s, r — > 0, 
while keeping the ratio s/r constant. We can now define the approximate solutions u s = 
u s (t, 6) to the Cauchy problem for the system (|3.3|) - (|3.6|) associated with the initial data 

u (d):=u(to,e), OeS 1 . 

It is convenient to avoid the discussion of boundary conditions by assuming that Uo is 
defined on the real line R, after periodically extending it beyond S . 

As usual, the scheme is defined inductively. First, the initial data are approximated by 
a piecewise constant function 

u s (t , 6) := u Q (6 h+1 ), 9 e [9 h , 9 h+2 ), h even. 

Then, if u s is known for all t < t^, we define u s at the level t = t\ as 

u s (t, 9) = u s (t-, 6 k/h+1 ), 9 e [9 h , 9 h+2 ), k + h even. 

Finally, the approximation u s is defined in each region 

< t < t\ + \, 9h-i < 9 < 9h+i, k + heven, 

from the approximate generalized Riemann problem 

^G\Us(tk/9h-i), u s {ti, dh+i); t k ,9hj, 

that is, 

u s (t,6) := u{t,9; u s {t k ,9} l -\), u s {t k ,9h+{); t k ,9nj, 

as introduced in (I5.2I). 
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5.2. Sup-norm and bounded variation estimates 

We will first establish that the sup-norm of the approximate solutions remains uni- 
formly bounded for all compact intervals [to, T] with £q > and T < +00 in the expanding 
case, and for all compact intervals [fo, T] with to < T < before the mass density blows up 
in the contracting case, that is, for all intervals [to, T] such that 

sup sup Jl s < C\(T) in the contracting case. (5.3) 

s [f ,T]xSi 

Subsequently, we establish a uniform bound on the total variation, again on compact time 
intervals. Recall that s and r must satisfy the standard stability condition 

- sup a s (£,•)< 1, (5.4) 
r s 1 

where we have observed that the function a s provides one with a bound on all characteristic 
speeds associated with the classical Riemann solver. In the expanding case, a s remains 
bounded by the supremum of the initial data flo but, in the contracting case, the bound 
imposes that the coefficient a s does not blow up. However, thanks to the equation (13.51 1, the 
bound i5.3} on the mass density implies an a priori bound on the function a s , so that the 
condition (|5.4|) can be ensured a priori from the initial data and once we know the constant 
d(T). 

Lemma 5.1 (Global sup-norm estimate). Consider any bounded initial data W°, V°,flo,/To,^o 
satisfying the physical constraints 

< fl , < Jl , \v \ < 1. 

Then, in the expanding case, for all T > to there exists a constant Co, depending on the time T, to 
and the initial data, and a constant C' Q depending on to and the initial data so that for all t e [to, T] 
and all relevant s 

0<l/Co(T)<fl s (£,-)<supflo(-), 

sup \W s + (t, •) - 1/(201 + 7 sup \V s + (t, -)| < C , 
s 1 1 s 1 



sup 

s 1 



( l + v s (t,-) \ 



+ 7-7W SU P I lo § ^(t, -)| < C (T). 

1+ K s i 



On the other hand, in the contracting case the same conclusion holds for any T e (fo, 0) such that 
d53b holds. 



Proof. We apply Propositions 14. 1 1 and 14.51 in which we derived "maximum principles" for 
the geometry and the fluid. Observe that it is at this juncture that it is important to formulate 
these principles in terms of the Riemann invariants associated with the geometry and the 
fluid variables. In each ODE step, the maximum of the geometric Riemann invariants Hf 
defined as in <(4.14)> does not increase in the expanding case, and may only increase by a 
factor 1 + C s in the contracting case. In each Riemann problem step, the variables H ± are 
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non-increasing as well. The fluid variables w, z are also similarly controled in both the 
ODE and the Riemann problem steps, while for the function a we need to require that 
the density is a priori bounded in the contracting case; that is, the assumption (|5.3|) . By 
iterating our estimates and covering a compact time interval [to, t], and after observing 
(1 + Cs) n < e-^-fo) provided (f - f ) < n s, we obtain the desired uniform bounds. 

In the contracting case, Proposition 14.51 shows that the amplification factor for Hf at 
each time step is 

t 2 k a s (t k+1 ,d) 2 

t 2 k+1 a s (t k ,e) 2 ' 

which we can control from the uniform bound (|5.3D , as follows, using (log a s ) t - -(1-k 2 ) t J[ s 
so that 

sup fls( ** +1 'f = sup exp ((1 - k 2 ) f M |f'| 6) df 

<exp^l-^)Ci(T)(tk + i-tid\io\j 

for all t k < tfr+i < T. By iterating this estimate over all times < T, we conclude that the 
sup-norm of Hf remains uniformly bounded in the compact interval [to, T]. □ 

We next establish the uniform bounds on the total variation of the approximate solu- 
tions. Denote by Uk,h+i the value achieved by the function u s at the mesh point (t^, dh+t), 
i.e. 

u k/h+1 := u s (t k , dji+i) k + h even, 
and denote by u k ^ the solution to the classical Riemann problem 

f^(Uk-l,hf u k-l,h+2'r h-lr 9h+l), 

where 

Jr>t k 
g(z, S T u k/h+1 )dT. 
k-i 

We divide the (t, 0)-plane into diamonds Ajy, (k + h even) with vertices (6k-i r h,tk-i): 
(6k,h-i'tk)> {9k,h+l' tk)/ (6k+i,h/tk+i)- To simplify the notation, we introduce the values of 
u s at the vertices of and the corresponding classical Riemann solvers by 

U S '■- u k-l,hf "W := W^/i-l/ U E '■- Wjyz+l/ U N '■- Uk+ifi, 

Uw '■= Uk,h-\> U E '■= Uk,h+l' U N := u k+l,h/ 
in terms of which, the strength e*(Ajy,) of the waves entering the diamond is defined as 

e*(A k/h ) := \e(u w ,u s )\ + \e(u s ,u E )\, 

whereas the strength e*(A^) of the waves leaving it is defined as 

e*(A Kh ) := \e(u w ,u N )\ + \e(u N ,u E )\. 
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Let / be a spacelike mesh curve, that is a polygonal curve connecting the vertices (O^h+i, h) 
of different diamonds, where k + h is even. We say that waves (wjt-i^, w^/2+1) cross the curve 
/ if / connects (6k-i,h' 6k-i) to (j/jyi+i/ h) an d similarly for ( Uk,h-lt Uk-x,h)- The total variation 
L(/) of / is defined as 

W) ■= ^ Huk-W^k,h+i)\ + \e(Vk,h-i,Uk-i,h)\, 

where the sum is taken over all the waves crossing /. Furthermore, we say that a curve / 2 
is an immediate successor of the curve }\ if they connect all the same vertices except for 
one and if J 2 lies in the future of ]\ . For the difference of their total variation we have the 
following result. 

Lemma 5.2 (Global total variation estimate). Let be two spacelike curves such that / 2 is 
an immediate successor of ]\ and. let Ajy, be the diamond limited by these two curves. Then, for each 
time T, strictly larger than the time ti determined by A^/or which the assumptions ofLemma \5.1\ 
hold, there exists a constant Co(T) depending on T and the sup bounds on the initial data, only, 
such that 

L(J 2 )-L(h)<C (T)se*(A Kh ). 

Proof. By definition, we have 

L(/ 2 ) - L(/i) = \e(u w ,u N )\ + \e(u N ,u E )\ - \e(u w ,u s )\ - \e(u s ,u E )\ 
= e*(^k,h) ~ £*(^k,h)- 

Observe that |e(ww/Wiv)l + |£(wn>we)I = \£(uw,ue)\ since i?iv is just one of the states in the 
solution of the Riemann problem for Uw,Ue- Hence, we can write 

L(J 2 )-L(J 1 ) = X 1 +X 2 , 

where 

X a := \e(u w ,u E )\ - \e(u w ,u s )\ - |e(M S ,u E )|, 
X 2 := \e(u w ,u E )\ - \e(u w ,u E )\. 
By Proposition l4.2l we have 

X 2 < 0. 

To estimate X 2 , note that, by Lemma |5?Ll u s is uniformly bounded for all t € [fo, T], hence 
by Lemma 1431 we can write for some constant C depending on T and the initial data 

X 2 < C \e(u w ,u E )\ (|ww - u w \ + \u E - u E \J + C |(w w - u E ) - (u w - u E )\ 

< Cs \e(u w ,u E )\ ( sup \u' w (t)\ + sup \u' E (t)\) + Cs \(u w - u E )\ 

te[t k -i,t k ] te[t k - u t k ] 

< Cs\e(u w ,u E )\ < Cs(\e(u w , u s )\ + \e(u s ,u E )\), 

since \(u w - u E )\ = 0(l)\e(u w ,u E )\. □ 
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5.3. Convergence and consistency properties 

In view of the uniform total variation bound, Helly's theorem allows us to extract a 
converging subsequence and, hence, to arrive at the following conclusion. 

Theorem 5.3 (Convergence and existence theory for the Einstein-Euler equations). Let 

uq - (W°, Vl,ao,JIo,vo) be an initial data set for the system d3.3IH3.6l) , which is assumed to 
belong to BV(S 1 ) and satisfy the physical constraints 

< do, < Jl , \v Q \ < 1. 

Consider the approximate solutions u s - (W s ± ,V s ± ,a s ,fi s ,v s ) constructed by the random choice 
scheme in Section 15.11 Then, the functions u s (t, •) belong to BV(S l )for each fixed time t e [to, +°o) 
in the expanding case and for each t e [to, T] in the contracting case, where T < is such that the 
mass-energy density satisfies the uniform bound 05.31 ) for some constant C\{T). Moreover, on any 
compact time interval, the solutions u s are Lipschitz continuous in time with values in L 1 , and have 
uniformly bounded total variation, that is, for all t, t' e [to, T] 

TV(u s (t,-))<Co(T)TV(u ), (5.5) 

f \u s (t, 0) - u s {t', 0)| d6 < C (T) TV(u ) (\t - f\ + s), (5.6) 

where Co is a constant depending on to, T and the supremum bounds on the initial data and, 
additionally, on Ci(T) in the contracting case. Furthermore, there exists a subsequence of u s 
converging in L 1 , whose limit u = (V, W,a,Ji,v) is a weak solution to the system ( |3.3b — (13.6b and 
satisfies the physical constraints 

< a(t, 0), < p(t, 6), \v(t, 0)| < 1 

and periodic boundary conditions. 

Proof. Lemma [5721 implies that we can estimate the total variation L(/^ +1 ), along the curve 
Jk+i connecting vertices of the form (0fc+i,/t, ifc+i) to vertices (0fc+2^+i/ tk+l)/ i n terms of the 
the total variation L(Jfc) of the curve fa. Hence, we have inductively 

L(fa +l )<L(fa) + sC(T)Y j e*(A k , h ) 

h 

< (1 + s C(T))L(fa)<e c ^L(h), 

and also 

L(h)<CTV(u ), 

which imply the claim made on u s . Observe that by periodically extending Uo to the 
real line, the proof of (|5.6|) proceeds along the same lines as the standard proof since, 
by construction, u s are periodic as well. By Helly's compactness theorem, there exists a 
subsequence, still denoted by u S/ converging strongly in the L 1 norm toward a function u 
with bounded variation. Moreover, the uniform estimates valid for u s extend to the limit 
u. 
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It remains to show that u is a weak solution to the system d3.3M3.6D , which we express 
in the form 

d f H(u, t) + d e F(u, t) + G(u, t) = 0, 
with obvious notation. By fixing a compactly supported smooth function (p, we have 

X+oo /->+oo 
I (d t (p H(u s , t) + d g <p F(u s , t) + <p G(u s , t))dOdt 

r(k+l)s Mh+\)r 

= Tj Tj \ (d t (pH(u s ,t) + d e( pF(u s ,t) + ci)G(u s ,tj)dOdt 

k h+kodd Jks J ( h ~ 1 '> r 



where 



r(k+l)s p(h+l)r 

Qs := Y" I I (pG(u s , t)dOdt, 

kh Jks J(h-l)r 

h+kodd 

[H(u s ((k + l)s- 0), (fc + l)s - ) <P((k + 1) S/ 9) 
(h-l)r \ 

- H(u s (ks+, 6), ks + ) (j)(ks, 0)J dd, 
[F(u s (t,(h + l)r-),t)(t)(t,(h + l)r)dt 

kh Jks \ 

- F(u s (t, (h - l)r+), t) cf)(t, (h - l)r) dt. 



k,h 
h+kodd 



and 



CT := 



k,h 
h+kodd 



For the first term, we have 

o^ott)^ L ( s2 + ? ' 2 )(s + '' + K,, I+ i-^- 1 i)^upp^ 

k h+k odd 

where ^supp^ denotes the characteristic function of the support of <p; see |Q] for further 
details. Since u s have uniformly bounded total variation and <p is smooth, we have Q} s 



for s — > 0. The second term can be rewritten as 

-(/i+l)r 



(H(u s (ks- %,)/ *s-) - H(w s (fo, ftr), fo))^(fo, 

-»+oo 



X+oo 
H(u s (t Q ,e),to)(t>(to,e)de / 
CO 



and it is well-known that the first part on the right-hand side goes to zero as s — > 
for almost every equidistributed sequence (fljt)iceiN- h remains to estimate Qg, but after 
rearranging some terms, we obtain 

X(k+l)s 
(F(u s (t, (h + l)r+), t) - F(u s (t, (h + l)r-), t))<f>(f, (h + l)r)df = 0, 
A' //+,voaa 1 
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since u s is smooth there. Collecting the previous estimates, we conclude that, as s — > 0, 

X+oo r*+oo 
J (d t (f) H{u s , t) + d e (pF(u s , t) + (j)G(u s/ t))dddt 

X+oo 
H(u s (t Q ,6),to)(p(t ,6)dd 
oa 

tends to zero. Hence, by the Lebesgue's dominated convergence theorem, passing to the 
limit u s — > u, we arrive at 

/->+0O /->+0O 

' J (<9 f ^ H(m, + d g cj)F(u, t) + 0G(u, tfjdddt 

X+oo 
H(u Q (t ,e),to)(p(to, 6>)d<9 = 0. 
oo 



□ 

6. Future expanding Einstein-Euler spacetimes 

We are now in a position to establish our first main result. 

Theorem 6.1 (Expanding Einstein-Euler spacetimes with Gowdy symmetry). Consider any 
BV regular, Gowdy symmetric, initial data set on T 3 for the Einstein-Euler equations, and assume 
that these initial data set has constant area to > and is everywhere expanding (toward the future). 
Then, there exists aBV regular, Gowdy symmetric spacetime M with metric g a p and matter fields 
p and u a satisfying the Einstein-Euler equations (|1.1[) - (|1.3[) in the distributional sense, and the 
following properties hold. This spacetime is (up to diffeomorphisms) a Gowdy-symmetric future 
development of the initial data set and is globally covered by a single chart of coordinates t and 
(6, x, y) e T 3 , with 

M = [t ,+ M )xr 2 , 

the time variable being chosen to coincide with the area of the surfaces of symmetry. 

This global existence result was established earlier by LeFloch and Rendall [16] for a 
different class of regularity, which provides less regularity properties on the spacetimes. 

Proof. Given an arbitrary initial data set u = (W^, V° ± , a^, pQ, v$), we have constructed (in Sec- 
tion |5j} a sequence of approximate solutions for the "essential" Einstein-Euler equations 
03.3ft — (13.6b and this sequence was found to converge to an exact solution u = ( W+, V+, a, JI, v) 
with BV regularity to the Einstein-Euler equations understood in the sense of distributions; 
Cf . Theorem 15.31 above. More precisely, u is locally Lipschitz continuous in time and has 
bounded variation with respect to the space variable, so for any fixed T € (to, +oo) and all 
t,t'e[t ,T] 

TV(u(t,-))<C (T)TV(u Q ), 

\u(t, 6) - u(f, 0)| de < C (T) [\t - t'\ + s) TV(u ). 

Once a solution to the first-order system (13.3lHl3.6l) has been found, we only have to recover 
the original metric coefficients U, A and v and show that the constraint equations remain 
satisfied for all times. However, this has precisely been done in Section 13.21 hence, the 
proof of Theorem l6.1l is now complete. □ 



X 
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7. Future contracting Einstein-Euler spacetimes 



7.2. Global geometry and energy junctionals 

We are now in a position to state our second main result. Observe that, for the con- 
tracting spacetimes under consideration now, a partial existence result was established in 
LeFloch and Rendall [16] for a different class of regularity, as it was already established 
therein that the development below covers [to, t c ) with t c < 0. 

Theorem 7.1 (Contracting Einstein-Euler spacetimes with Gowdy symmetry). Consider 
any BV regular, Gowdy symmetric, initial data set on T 3 for the Einstein-Euler equations, which 
has constant area -to > and is everywhere contracting toward the future. Then, there exists a BV 
regular, Gowdy symmetric spacetime M with metric g a p and matter fields \x and u a satisfying the 
Einstein-Euler equations l|l.l|) - (|1.3|) in the distributional sense, and the following properties hold. 
This spacetime is (up to dijfeomorphisms) a Gowdy-symmetric future development of the initial data 
set and is globally covered by a single chart of coordinates t and (6, x, y) e T 3 , with 



for some t c < 0, the time variable being chosen to coincide with minus the area of the surfaces of 
symmetry and so that the mass energy density blows-up at t c , that is, 



the future development exists up to the maximal time t c = when the area of the surfaces of 
symmetry shrinks to zero. 

The theorem above seems to be essentially optimal and, for instance, we do not expect 
that a global-in-space, energy-type functional should allow us to decompose the space of 
initial data into spacetimes that blow-up (that is, tc < 0) and spacetimes that do not blow- 
up (that is, tc - 0). At this juncture, let us recall that the class of T 2 -symmetric vacuum 
spacetimes was studied by Isenberg and Weaver who could exhibit an energy that 
vanishes precisely for a set of exceptional solutions and also controls the behavior of the 
area function. Later, Smulevici I28T) was able to handle spacetimes satisfying the Einstein 
equation with a cosmological constant. In contrast to these problems, the Euler equations 
are nonlinear partial differential equations of hyperbolic type, so that the analysis of the 
global geometry of Einstein-Euler spacetimes is technically much more involved than 
the one of Einstein-vacuum or Einstein-Vlasov spacetimes. By taking advantage of the 
property of finite speed of propagation satisfied by the Euler equations, it is straightforward 
to "localize" our conclusions above and to construct solutions to the Einstein-Euler that 
are small perturbations in the L 1 norm (but not in the BV semi-norm) of non-blow-up 
solutions (for instance homogeneous spacetimes) and, however, blow-up in a time tc < 
with tc arbitrarily close to the initial time to. 



M = [t Q ,t c )xT 2 




(7.1) 



Furthermore, when the initial data set satisfies 




(7.2) 
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The rest of this section is devoted to the proof of Theorem 17. 1 1 and will make use of the 
energy functionals 

Ei(t) := J h dd, E 2 (t) := J (h + hf) dO, 



with 



(tW + - 1/2) 2 (fW_ - 1/2) 2 _^ M ft l + feV 

1 2«i 2 2flf 2 8«£ 2 8a£ 2 ' 1- a 1-y 2 ' 



Observe that /?i can be obtained from the standard energy density by a simple transforma- 
tion of the metric coefficient U \-t 211 — log t . Introduce the flux 



(JW+-1/2) 2 (£W_-l/2) 2 V% V z _ 



2af 2 2at 2 Sat 1 Sat 2 ' 

and observe that hi, g\ satisfy the balance laws (cf . with (|4.13|l ) 

d t hi + d e {agi) = —hi - 
dtgi + d e (ahi) = jgi - y> 

withfci := ^ W++ ^ t 2^ ^ + ^160^ • By using the evolution equations, we obtain the following 
result. 

Lemma 7.2. Both functionals E\ and E 2 are monotonically increasing and, specifically, 



§«> = 



dt ^ t f s i 



(f(W + + W_) - l) 2 (V + + V.) 2 M ^3k 2 + l 



2flf 2 Sat 1 



+ — — 7^+hf + --— — \dd>0. (7.4) 



Hence, we have < —^r- and thus 



dt 

Ei(t) < E 2 (t) < E 2 (t )j, 

which implies that the functionals E\ and E 2 do not blow up before the singularity hyper- 
surface t = is reached. 

Observe that the arguments in Section [6] can be repeated, once we establish that the 
metric coefficient a and fluid density y. do not blow-up. The new difficulty comes from 
the uniform bound (|5.3|) on the mass-energy density. In fact, as discussed now, there are 
solution where this condition does not hold uniformly up to t - 0. 
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7.2. Revisiting the class of spatially independent solutions 

Before we provide a rigorous proof of our non-blow-up result, we provide a simple 
derivation for special classes of solutions. By inspection of the proof of Lemma 14.41 one 
easily checks the following slightly stronger statement. 

Lemma 7.3 (No blow-up property for the homogeneous solutions in the contracting case). 

Fix a time fo < and initial data Uq = (Wq , V*, uq, Jlq, vq) satisfying the constraints 

0<a , < p , \v Q \ < 1, (7.5) 



and consider the corresponding solution u = u(t) of the ordinary differential system (14.51) , defined 
on its maximal interval of existence [to, T) with T < 0. Then, for the geometric variables, one has 



W ± (t)\ < C ^P-, \V ± (t)\ < Q >/«(*), «o < a(t), 



while, for the fluid variables, 

ll + v 

log 



2k a(t) 
llogfil < Co—, 



1 + k 2 ' 

where the constant Co is uniform in time and depends on the initial data to, uq, only. 

The above lemma provides a quantitative control of all the variables (but a) in terms of 
the coefficient a. If a(t) blows-up when t approaches some critical time T e (to, 0), then at 
that time the density may tend to or +oo. To be more precise and describe the nature of 
the blow-up for the fluid variables, we can use the energy 

E _ (tW+ - 1/2) 2 + (fW_ - 1/2) 2 + v l + v - 



2at 2 2at 2 8at 2 8at 2 ' 

introduced in the previous section, as a criterion for the blow-up of a. 

Lemma 7.4. Given any initial data uo and the corresponding solution u = (W ± , V ± , a, Jl, v) of the 
ordinary differential system d4.5l >, one can distinguish between the following two cases: 



IfEi(to) + 0, then the function a remains bounded. 



• IfE\(to) - 0, then the velocity remains bounded away from ±1 and the following two subcases 
arise: 

- When fio is sufficiently large and, specifically, 

3 f o Mo > 1/ 

the coefficient a and consequently the density fi ^ blow-up for some T € (to, T). 
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- When /io * s sufficiently small and, specifically, 

then a, Jl remain bounded. 

Our latter claim is just a special case of the final statement in Theorem 17.11 Observe 
that, in the spatially homogeneous case, where, in addition, v = 0, the system can be solved 
explicitly and the above lemma actually covers all cases. 

Proof. If Ei (to) > 0, we have 

(a-'EiY = -2^3 + W_) - l) 2 + (V + + y_) 2 /4) > 

and, therefore, a~ l (t) > a' 1 Ei(to)/Ei(t), so that the coefficient a is bounded. On the other 
hand, if the energy initially vanishes Ei(io) = 0/ by uniqueness, it must remain identically 
zero for all times. Therefore, the system (|4.5I) simplifies and 

tW = 1/2, V = 0, a' = (1 -fc 2 )!^^, 

with W := W+ = W_, V := 7+ = V-, and 

^-l+AVV ,9,3/c 2 + l 



/ , _1 +]fcV\ i l2 ' 



4f 



! -(1 + k 2 )v 
After changing the variables, we obtain 



/, (3/c 2 + l) l + i> 2 

(logfl- 1 !^) - 



and, by integration and using \v\ < 1, 



4|f | (1 - k 2 v 2 ) 



. ,(3fc 2 +l)/4 

to 



a 1 tjj< a lf oMoly 

or 

4 // f f 1 "^ 4 



Consequently, our assumption fij^o > 1 implies that there exists a critical time t e (to,0) 

a j~ / /t \3(1-^)/4n 

such that 1 = |fp^o(l _ ( L which leads us to a (t) = 0, a contradiction. The latter 
claim in the lemma can be checked similarly. □ 
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7.3. Estimates on the geometry and fluid variables 

We now consider any solution to the Einstein-Euler system defined on some interval 
of [to, t c ) and we derive a bound on sup s i a which depends upon a lower bound on t h 
J sl a~ 1 (t,6)d6. This property will provides us with the key ingredient of the proof of 
Theorem l7.ll 

Proposition 7.5. If there exists a function C = C(t) defined for t e [to, t c ) such that for 11 relevant 
times t 

0<C(t)< | a- 1 (t,0)dO, (7.6) 

then there exists also a function C\ - C\{t) such that for all such times 

•»t 77 



sup f -H—fr^zZd®. (7.7) 

M 

Then, a fortiori, sup s i J |u(t, •) dx is bounded and the metric coefficient a is uniformly bounded on 



To establish this result, we rely on the following two lemmas, which we derive first 
and concern the infimum of the time-average of the mass density ^ 777^2(1", ~)dz and its 
variation in space, respectively. 

Lemma 7.6 (Lower bound for the averaged mass-energy density). For all t e [to, t c ), one has 

(1 + k 2 ) I a~\t, 6)de ] j inf J JL ( T , .) dT < |t | (E 2 (t) - E 2 (t )). 
Proof. Namely, from the energy-type estimate (|7.4[) we deduce 

X X/'^TT^r^ ^ ^ * l*ol (£ 2 (0- E 2 (t )). 

Since fl -1 is monotonically decreasing in time, we obtain 
f a-^^e^inf f -l—(r,.)dx 

< r a- 1 ^^) r ~t!—(T,9)dTde< r r ^JL^dBdz. 

JS 1 Jt 1 - " Jfn Jsi 1 - " 



□ 



Lemma 7.7 (Variation in space of the averaged mass-energy density). For all t e [to, t c ) and 
60, Q\ £ S 1 , one has 

J to Jfo 1 - V 1 

+ |f |(Ei(t) - Ei(f )) + 2\t\E 2 (t) + 2\t \ E 2 (t ). 
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Proof. Setting M := > 0/ integrating the "second" Euler equation in (I3.6D over some 
slab [to,t] X [0O/6i]/ an d finally using |u| < 1, we obtain (recalling that the time variable 
takes negative values) 

(l + /c 2 ) f xM^eOdT-fc 2 f tM(t,0 o )^t 
< I tM (fc 2 + ^ 2 )(t, 0\)dx — I zMik 2 + v 2 )(T,6 )dT 

J tg Jtg 

t— g! dOdz + (1 + fc 2 ) (te _1 Afo(f, 0) - toa^Mvito, 9)) d6. 
a a Je 

In the above identity, the double integral term is controlled by E\, i.e., 

f f 1 ^z gl d9dz < - f f ^xhxdddx < |f | (Ei(f) - Ejfo)), 

whereas, for the other two terms of the right-hand side, 

(l + /c 2 ) (ta _1 Mu(t, 0) - toa^Mvito, 9)) d6 
< t 2 E' 2 (t) + t 2 E' 2 (t ) < -2t E 2 (0 - 2t E 2 (t ). 

Hence, we obtain 

-k 2 f tM(t, 9\)dz 
J t 

<-(l + k 2 ) [ TM(T,9 )dT-to(E 1 (t)-E 1 (to))-2tE 2 (t)-2t E 2 (to), 

J t 

which is the desired estimate. □ 
Proof of Proposition 17.51 Combining Lemmas |7.6| and [7!7l we obtain 

^sup f B^(r,.)dr<(l + ^)ird f I^-(r r )dz 
# Jto l-v 2 si Jt 1-v 2 

\t \ (Ej(f) - EM) + 2\t\ E 2 (t) + 2\t \E 2 (t ) 
C (t)^j^ a- 1 (t,9)d9^j +Ci(f) 



•»t 771 
si Jt 

+ 

< 



for some (explicit) functions Co(t), Ci(t) > 0. □ 

Proof of Theorem [77TI We now complete the proof of our second main result. In view of 
Proposition |7.5| and the earlier results in Section 4, it is sufficient to derive the lower bound 
(|7.6[) on J^j a -1 dd. Specifically, we now check that if the time interval [to, t c ) is sufficiently 
small or if the initial data set satisfies the condition (|7.2|) then (for all t e [to, t c ) or else for 
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all t £ [to, 0), respectively) there exists a function C(t) > such that 07.61 ) holds. Namely, by 
integrating the "first" Euler equation in d3.6l > over [fo, t] X S 1 we obtain 

r r _ 1 + k 2 v 2 

Js 1 Js 1 1 - ^ 

r , 1 + fc 2 y 2 r ( r ., 3A- 2 + 1 
< «- %lf |- ^d9 + a~ l Ji——dedT 

and, by Gronwall's inequality, 



r it \ {3k2+m 



with 



C(t )= f AoVolfoKi+^A 1 -^)^- 

Js 1 

Therefore, in view of (fl _1 )f = a _1 l ui(l - k 2 ), we arrive at 

£ fl- 1 ^ > ^ fl-He)^ + 1 c(fo)|t | ((^) 3(1 "" 2)/4 

which implies (I7.6I) as long as the right-hand side above remains positive, that is, either on 
a sufficiently small interval [to, t c ) (since the lower bound is always positive if t is close to 
to), or else on the whole interval [to, 0) if the small mass condition (|7.2[) holds. We conclude 
by applying Proposition |7.5[ which provides us with the required uniform bound in order 
to complete the proof of Theorem 17. II □ 
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